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ABSTRACT 

Measurements of cosmic microwave background (CMB) anisotropics by interferometers offer several 
advantages over single-dish observations. The formalism for analyzing interferometer CMB data is 
well developed in the flat-sky approximation, valid for small fields of view. As the area of sky is 
increased to obtain finer spectral resolution, this approximation needs to be relaxed. We extend 
the formalism for CMB interferometry, including both temperature and polarization, to mosaics of 
observations covering arbitrarily large areas of the sky, with each individual pointing lying within 
the flat-sky approximation. We present a method for computing the correlation between visibilities 
with arbitrary pointing centers and baselines and illustrate the effects of sky curvature on the £-space 
resolution that can be obtained from a mosaic. 

Subject headings: cosmic microwave background — techniques: interferometric 



1. INTRODUCTION 

The study of anisotropies in the cosmic microwave 
background (CMB) radiation has revolutionized cosmol- 
ogy. Key to this revolution have been coupled advances 
in theory, data analysis, and instrumentation. In partic- 
ular, the design of experiments with exquisite systematic 
error control has been crucial for progress in the field. 
Interferometers offer several advantages in this respect, 
with simple optics, instantaneous differencing of sky sig- 
nals without scanning and no differencing of detectors. 
The shape of the beam can be well understood and the 
measurement is done directly in Fourier space where the 
theory most naturally lives. 

Pioneering attempts to detect C MB anisotropy with 
interfe romete rs were made by iMartin fc Partridge! 
(|1988J) and ISubrahmanvan et alJ l|1993j) . Several 
groups have successfully det ec ted p ri mary CMB 
anisotropies iO'Sullivan et all 119951: iBaker et al l 
19991 IHalverson et all I2002t IPearson et all 12003 
Tavk3r~et^uTT20"03l) and polarization jReadT^adet^Ll 
20041: iLeitch et all 12005(1 . using interferometers. The 
formalism for analyzing CMB data from int e rferom - 
eters has been devel oped by IHobson et alJ 1 1995 ): 
Hobson fc Magueiiol lll996tk iWhite et all 1 1999): 
Hobson fe Maisingerl l)2002ft : and iMvers et alJ (|2003|i ; 
as well as in th e expe rimental pa p ers c ited above. 
iPark et all <H)03) and iPark fc Nd l|2004fl examined 
interferometric polarimetry. 

In the Fraunhofer limit an interferometer measures the 
Fourier transform of the sky, multiplied by the primary 
beam. The primary beam determines the instantaneous 
field of view of the instrument and its Fourier transform 
is simply the autocorrelation of the Fourier transform 
of the point response of the receiver to an electric field. 
The angular scale probed by any pair of telescopes being 
correlated is determined by their spacing in units of the 
observational wavelength 1 . The range of scales probed 

1 We assume throughout monochromatic radiation; the general- 



by the interferometer is then determined by the spacing 
of the elements, while the resolution in spatial wavenum- 
ber is determined by the area of sky surveyed. By "mo- 
saicking" several smaller patches together, the resolution 
in spatial wavenumber can be increased, although the 
range of spatial scales remains fixed by the geometry of 
the interferometer elements. 

In most cases it has been assumed that the field of 
view is small, so that one can use the "small-angle" 
or "flat-sky" approximation. However, if we want fine 
resolution in spatial wavenumber - which future exper- 
iments are d riving towards - we need to survey large 
areas of sky l|Hobson fc Magueiiol 1 1996j) and thus relax 
this assumption. The purpose of this paper is to extend 
the formalism presented in the above papers to the case 
where each individual pointing of the interferometer is 
still within the flat-sky approximation but by mosaick- 
ing many pointings together a significant area of sky is 
surveyed. Our extension allows one to see how large an 
error is being made in assuming the flat-sky approxima- 
tion and shows how corrections can be systematically 
incorporated. 

The central idea of this paper is the following. The 
key ingredient in analyzing a mosaic of interferometer 
pointings is the set of two-point visibility correlations. 
For each pair of pointings, we can calculate the correla- 
tions in a spherical coordinate system that places both 
pointing centers on the equator. If each pointing has a 
small field of view, then we can approximate the sphere 
by a cylinder in the vicinity of the equator, allowing the 
use of Fourier analysis rather than a more cumbersome 
expansion in spherical harmonics. 

The outline of this paper is as follows. We begin in fJ5] 
by reminding the reader of some basic results in the flat- 
sky limit. We then show how this can be extended using 
a cylindrical projection in J^and make contact with the 
exact spherical harmonic treatment in <2J Section [5] cx- 

ization to a specified frequency band is straightforward. 
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tends our results to include polarization, and we conclude 

in jn 

2. FLAT-SKY LIMIT 

We begin by considering the flat-sky limit and focusing 
on temperature anisotropics. Thus our interferometer is 
measuring a scalar field, T(x), defined on the 2D plane. 
In this limit the fundamental observable, a visibility, can 
be written 



^(u)~ Jd 2 x AT(xM(x)e 2 



(1) 



where dB v /dT converts from temperature to intensity 
units and A(x) is the primary beam (typically normalized 
to unity at peak) . From now on we will neglect the flux- 
temperature conversion factors and write T for AT. 

For Gaussian fluctuations, such as the primary CMB 
anisotropies, we need to compute the visibility correla- 
tion matrix 



V; 



= (V^V*^)) , 

where Ui and represent the baselines to be correlated 
and (•) represents an ensemble average. This can be re- 
lated to the usual temperature correlation function 

(r(m)T(n,-)) = r £(2* + ■ nj) (3) 

1=2 

for temperatures measured in directions and r\j where 
Cg are the multipol e moments. In our flat-sky limit, for 
a sing le patch fe.g- IWhite et al.l H999) 



(2) 



Vij cx j d 2 wS(w)A*{2n[w + u t })A(27r[w + u 3 }), (4) 

where the angular power spectrum S(u) is defined by 

(2tt) 2 u 2 S{u) ~ 1(1 + l)d for £ = 2nu, (5) 
and A is the Fourier transform of the antenna pattern, 

I(k) = (2vr)- 2 / d 2 xA(x)e- 1 ^. (6) 



The extension to multiple different patches, each with its 
own pointing center V, merely inserts a phase factor 

V^ lP2 cx J d 2 «;S , Hl*(27r[w+u J ])l(27r[w+u J ])e 27riW - D , 

(7) 

where D is the separation between the pointing centers 
V\ and Vi ■ From now on we will drop the superscripts 
on Vij . 

3. CYLINDRICAL METHOD 

The flat-sky approximation above is valid only if both 
the field of view of an individual pointing and the separa- 
tion D between pointings are small. We will now assume 
a mosaic of pointings that cover a large area, although 
each individual pointing observes only a small patch of 
sky. We therefore relax the second assumption while re- 
taining the first. 

For a statistically isotropic temperature field we are 
free to use any coordinate system we like to compute the 
visibility correlation V12. In particular, we can arrange 
to have the two pointing centers lie on the equator of a 
spherical coordinate system (8 — 7r/2) and be separated 
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Fig. 1. — Window functions for the covariance between two in- 
terferometer pointings in both the fiat and cylindrical approxima- 
tions. The antenna pattern is Gaussian with beam width a = 5°. 
The two pointings are separated by 120°. Each visibility has a 
baseline of magnitude u = 20 pointing in the </> direction. The 
plus signs are the cylindrical approximation, and the stars are the 
flat approximation. For clarity, only the real parts of the window 
functions are shown. 



by an angle f3. We introduce a cylindrical coordinate 
system, with the cylinder tangent to the sphere at the 
equator, denoted by £ = ((f>,z). Since both observations 
sample only regions near the equator, we can pretend 
that the data live on the cylinder rather than on the 
sphere. In this approximation, it is natural to expand 
the temperature T(x) in a discrete Fourier series in <fi 
and a continuous Fourier transform in z: 



T(i) = J2j dn 7 1 m(w)e l(m0+ ™ 2:) 



(8) 



with 



T m (n)T m ,{n ) J = — ^ 8 mm >o(n -n). (9) 

The power spectrum is V(k) = (2ir) 2 S (k / 2ir) with S as 
defined in Eq. J5j. It is related to the spherical harmonic 
angular power spectrum by V(k) ~ Ck for large k. 
The visibility becomes 



V(Vx,m) - (2tt) 2 J dnJ2T m (n)A* (2tt[ui + w]) . 



where the vector w has coordinates (w<p,w z ) 
(m, n) / (2ir) and 



(10) 



A(k) 



d 2 a 

(2ny 



-ik-£ 



(11) 



is the usual Fourier transform of the primary beam. Since 
we imagine A is non-zero only over a small region we can 
extend the integral over the entire plane. 

The visibility for the second pointing center is analo- 
gous, except for a phase factor e lm/3 , so the correlation 
between two pointings becomes 



V12 oc J dn^2 S , (w)yVi2(ui,u 2 ,w)f 



im{3 



(12) 
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Fig. 2. — Pointing centers (+) and baseline vectors (horizontal 
bars) for the mosaicking example described at the end of section 
Fil In calculating the covariance between the two visibilities at the 
upper corners of the grid, we must rotate to a coordinate system in 
which the great circle connecting them (the dashed line) becomes 
the equator. Note that the two baseline vectors are not parallel in 
the new coordinate system. 



Fig. 3. — Improvement in £-space resolution due to mosaicking. 
The solid curve is the window function for a single pointing. The 
dashed curve is the window function for the sum of all pointings 
in a 10 X 10 grid, after correctly accounting for baseline rotation. 
The dotted curve is the window function that would be obtained 
by incorrectly assuming a flat sky over the entire mosaic. All three 
window functions have arbitrary normalization. 



where we have defined the window function 

W 12 ( Ul ,u 2 ,w) = l*(27r[u 1 +w])2(2 7 r[u 2 +w]) (13) 

It is convenient to define a window function that is aver- 
aged over direction: 

V 12 = ^T*f 2) Q. (14) 

t 

( 12") 

To compute W t , we divide the integral and sum in 
equation (fT^fl into bands with £ — \ < 2nw < I + h. 
Within each band we assume the power spectrum is con- 
stant and remove it from the integral. We can calculate 
this window function in the flat-sky approximation in- 
stead of the cylindrical approximation if we like, simply 
by replacing the sum over m by an integral. 

Figure n illustrates the difference between the flat and 
cylindrical approximations for large (3. The difference is 
most significant when the baseline vectors Uj are nearly 
equal in magnitude and parallel to the separation di- 
rection (/>: otherwise the correlation V12 is always small 
whenever /3 is large. 

We can use this prescription to calculate the full visi- 
bility covariance matrix for a mosaic of many pointings. 
For each pair of pointings, we must transform to a coor- 
dinate system in which both pointings lie on the equa- 
tor. In performing this rotation, the components of the 
baseline vectors Ui will naturally be transformed. The 
Appendix contains an explicit recipe for performing this 
rotation. 

Figures and 01 present a simple illustration of how 
mosaicking increases the ^-space resolution of an experi- 
ment. In each pointing, the beam pattern is a Gaussian 
with beam width a — 5°. We assume a 10 x 10 grid of 
pointings, separated by 5° in both 9 and in a spheri- 
cal coordinate system, with the center of the grid on the 
equator (9 = 7r/2). For each pointing, we consider only 



a single baseline u = 220. Figure [21 shows the locations 
of the pointing centers and the baselines in Aitoff projec- 
tion. Note that although all the baselines have identical 
components in the spherical coordinate system, they do 
not when rotated to the appropriate coordinate system 
for computing the covariances. As an example, to com- 
pute the covariance between the two pointings in the up- 
per corners of the grid, we must use a coordinate system 
in which the great circle represented by the dashed line 
becomes the equator. In this coordinate system, the two 
baseline vectors have 6 components of opposite sign. 

The visibility obtained from any single pointing pro- 
vides an estimate of the power spectrum with a fairly 
wide window function (solid curve in Fig. |5J. We can 
obtain an estimate of the power spectrum with a nar- 
rower window function by simply adding together all 100 
visibilities. To find the window function for the sum, we 
write down the absolute square of the sum of all 100 
visibilities: 

/ 100 2 \ 100 100 

(E^ )=EW> = E V «- ( 15 ) 

\ i— 1 / 1 1 

The window function for the sum of all the visibilities 
is therefore the sum of w[ over all visibility pairs i, j. 

(a) 

We can compute each W z using the recipe described 
above. Specifically, for each pair of pointings, we find a 
rotation that brings both pointing centers to the equator, 
apply that rotation to the vectors Ui , u 2 using the recipe 
in the Appendix, and apply equation l|13() . 

The dashed curve in Figure is the window function 
for the sum of all pointings. As expected, the £-space 
resolution has improved. The dotted curve shows the 
window function obtained by incorrectly assuming the 
sky is flat over the entire mosaic - that is, assuming that 
all of the baseline vectors illustrated in Fig. [5] lie in the 
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same plane and are parallel. 

The difference between the dashed and dotted curves 
is almost entirely due to "baseline rotation" - the fact 
that, e.g., the baseline vectors at the corners of Fig. 
are not in fact parallel to each other when viewed in 
a coordinate system in which both lie on the equator. 
It makes virtually no difference whether we use the flat 
method (integral over to), the cylindrical method (sum 
over to), or an exact spherical harmonic calculation as 

described below in calculating each W* , as long as we 
get the baseline rotation right. 

Of course, other linear combinations could be used in- 
stead of a simple sum of all 100 pointings, resulting in 
window functions with peaks in different places (within 
the envelope set by the single-pointing window function). 

4. HARMONIC METHOD 

We can also make direct contact with usual spherical 
harmonic treatment in which 



and 



The visibility for a single pointing is 



(16) 
(17) 



V(u) = I d 2 rA(r)T(r)e 2 ^ = ^a /m F /m (u), (18) 

i.m 



where 



F im (u) = / A(v)Y lm (i)e 



27T2U-r 



(19) 



It is of course possible to perform these integrals numeri- 
cally and calculate the visibility covariance matrix with- 
out any approximations at all. In this section we will see 
how to obtain the cylindrical approximation from this 
exact expression. 

In previous treatments (e.g.. IWhite et all fl9991 . the 
flat-sky limit for a single pointing was taken by approx- 
imating the spherical harmonics near the pole of the 
spherical coordinate system (9 = 0). To obtain the vis- 
ibility covariance for two different pointings it is more 
convenient to place the pointing centers on the equator 
as in the previous section. 

Near the equator z = cos 8 ~ ir/2 — 9. Using the re- 
curr ence relations for the associated Legendre polynomi- 
als lAbramowitz fc Stegunlll972t iGradshtevn fc Rvzhikl 
19801), one can show 



Ye m (cf>,z) 
with 



Nt m e 



net, 



im.4, / cos ne m z if I - 
-i sinri£ m z if t ■ 



= y/i{t + 1) - m 2 



to even 
to odd 



(20) 
(21) 



and 



l]/2)! 



(i + m)\ {[t-m]/2)\ 



(22) 



which can also be written 



^ m =(-i) ( ^ +m)/2 2-y 



2£ + 1 

4"7T 



y/(e + m)\(e-m)\ 
by using ^radshtevn fc RvKhildlTflS^ 



»+-.! 



(2n + l)! 
2 2n+i n! ■ 



(23) 



(24) 



In all of these expressions the factorials should be inter- 
preted as r functions for non-integer arguments. In the 
limit when all of the fa ctorial moments are large we can 
use the approximation l)Abramowitz fc Stegunlll972l 

In Nl ~ (n + M In N - N + const (25) 

to write the normalization factor as 

(_ 1 )(I+m)/2 + \ 



Nir, 



TT (/2_ m 2)l/4 



(26) 



Note that the Yg m are eigenfunctions of the 2D Lapla- 
cian with eigenvalues — 1(1 +1), and the form of ni m 
guarantees that this is preserved in the cylindrical coor- 
dinates: 



d 2 



i:) 2 



Y f 



(f> 2 dz 2 
= -l{l+\)Y tm , 



I III 



(27) 

(28) 
(29) 

where the arrow indicates the cylindrical coordinate 
limit. 
If we define 



aim = (-1) 



l+m 



air, 



■aim, (30) 



then we can rewrite Eq. I|16(l as 

T(0 = J2^ m e l ^ + a em e 1 ^ 



(31) 



with the definitions k = (m,n^ m ) and k = (to, — ng m ). 
From now on we will take the sum over both ±n^ m as 
implicit and write 



T(£) = 5> ke ik '«. 



(32) 



This way of writing the spherical harmonic expansion 
makes the correspondence with the Fourier representa- 
tion explicit. 
The visibility becomes 

V(u) = (2tt) 2 ct(k)A* (k + 2ttu) , (33) 



and the correlation matrix is 
(2^) 4 



Vis 



^|A^ m |^>Vi 2 (ui,u 2 ,k/27r)e ik A (34) 



plus oscillatory terms that average to zero in the sum 
over I and to. 
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Fig. 4. — The Fourier modes = (m, ±n; m ). 



If we work at large values of £ and m we can replace the 
sum over k with ^ m and J dn and recover our cylindrical 
result, Eq. (|12|) . One can verify by explicit computation 
that using 

d£ _ yJt{l+\)-m\ 



dn 



£ + \ 



to turn the sum over £ into an integral and using the 
asymptotic form of Ng m in Eq. I|34|l leads to Eq. (|12|l . 

As Fig. 0] illustrates, the differences between the har- 
monic and Fourier expansions lie in the non-uniform grid- 
ding of the wave vectors k. For any given wavcnum- 
ber £ the modes are packed more sparsely near the m 
axis, corresponding to \m\ ~ £, and more densely near 
m ~ 0. The normalization factor Ni m weights the higher 
to modes more strongly to make up for this. 

Surprisingly, in some cases the cylindrical approxima- 
tion proves numerically more accurate than the approx- 
imation in this section, as illustrated in Fig. [5J Specif- 
ically, in cases where the baseline vectors point in the 
4> direction, the sum in equation (|34|) is dominated by 
modes with \m\ = £, which are precisely the modes for 
which the plane- wave approximation to Yi m is worst. We 
found no instance in which the cylindrical approxima- 
tion l|20(l does worse than the approximate spherical har- 
monic expansion of this section, so for numerical work 
one should either use the cylindrical approximation or 
the exact full-sky expressions. In general, we find that 
the cylindrical approximation starts to become poor for 
Gaussian beam widths cr ~ 8° (FWHM ~ 20°). 

5. POLARIZATION 

In this section, we will generalize the results of the 
previous sections to observations of linear polarization. 
Instead of considering a single scalar observable T, we 
must consider observations of the two Stokes parameters 
Q,U. 



It is convenient to combine the Stokes parameters into 
the complex quantities 



P+ 



1 



-j={Q±iU), 



(36) 



because these quantities transform in a relatively simple 
way under rotations: under a rotation by an angle "0 
about a given point f, P±(r) — > P±(r)e T2 ' 1 ^ . In other 
words, P + is a quantity of spin weight —2 and P_ has 
spin weight +2. These transformation properties make 
(P_l_ , P_ ) a more convenient basis of observables to work 
with than {Q,U). The two bases are related by a unitary 
transformation, 



Wi i 

a/2 V~* 1 



u 



(37) 



so any results derived in one basis can easily be trans- 
formed to the other. 

An interferometer that works by interfering circularly 
polarized radiation from the two antennas measures the 
visibilities 



V±(u) = / d 2 r P±(r)A(r)e 



27rm-r 



(38) 



Specifically, interfering left-circularly polarized radiation 
from antenna 1 with right-circularly polarized radiation 
from antenna 2 yields V+, and reversing the senses of 
both circular polarizations yields V_. (Interfering right 
with right and left with left yields visibilities that probe 
total intensity and circular polarization.) On the other 
hand, an interferometer that works by combining linear 
polarization states would measure visibilities Vq and Vjj 
for the individual Stokes parameters. (For instance, in- 
terfering E x from antenna 1 with E y from antenna 2 
yields Vu-) We will assume that the measured 2 quanti- 
ties are V± rather than Vq^Ui but all results are easily 
transformed to the Q, U basis using Eq. (|37|l . 

As in the case of temperature anisotropy, the key in- 
gredient in analyzing CMB interferometric observations 
of polarization is the visibility covariance matrix: 



V 



12 



(V r ±(ui,7'i)7±(ua,^i) 



(39) 



where and Pi represent baselines and pointing cen- 
ters for a pair of visibilities. Note that in this equation 
the signs of ± and ± can be varied independently - 
that is, there are in general four distinct covariances, 

v ++ v H — v — h 

^12 ' ^12^ i ^12 _J ^12 ■ 

Our primary interest will continue to be the case where 
the flat-sky approximation is appropriate for each indi- 
vidual pointing but the separation between pointings is 
not necessarily small. We will present exact expressions 
for the visibility covariances in terms of spherical har- 
monics first, then show that they reduce in the this limit 
to cylindrical-sky expressions similar to the anisotropy 
results above. 

5.1. Flat sky 

2 We are ignoring some sources of systematic error in this ex- 
pression. For instance, in an instrument with cross-polar beam re- 
sponse, each measured visibility would contain contributions from 
both and P_, with different effective antenna patterns. 
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Fig. 5. — Comparison of approximations. The window function 
is shown for the covariance between visibilities with u = 20(p, beam 
width a = 5°, and pointing centers separated by /3 = 15°. The 
solid curve is the exact spherical harmonic calculation; the dashed 
curve is the cylindrical approximation; and the dotted curve is the 
approximation obtained by approximating the spherical harmonics 
by plane waves. 



Assume our observations cover a small enough patch 
of sky that we can replace spherical harmonic expansions 
with Fourier transforms: 



P±(x) = / d 2 fcP±(k)e 



ik-x 



(40) 

The two Fourier transforms are related like this: 

Pi(k)=P T (-k). (41) 

A key insight into the analysis of CMB po l ariza- 
tion was the observation l)Kamionkowski et alJ Il997t 
IZaldarriaea fc Seliakl Il997|) that any polarization field 
can be decomposed into a scalar part (conventionally de- 
noted E) and a pseudoscalar part (denoted B). The E-B 
separation is particularly simple in Fourier space: modes 
with polarization parallel or perpendicular to k are E 
modes, while modes polarized at 45° are B modes. In 
terms of P±, this means that 

P±(k) = (£(k) ± iB(k))e ±2 ^ k , (42) 

where V'k is the angle made by the wavevector k with 
respect to the x axis. 

Assuming that the polarization was generated by a ho- 
mogeneous, isotropic, parity-respecting process, the two- 
point correlations between E and B are determined by 
two power spectra Ve,b 

(E(k)E*(k')) = (2ir)- 2 V E (k)S(k - k'), 

(P(k)P*(k')> = (2ny 2 V B (k)6(k - k'), 

(£(k)P*(k'))=0. 

This means that the covariances of the polarization P± 
are 

(P± (k)Pl (k')) = (2tt) - 2 (T E (k) + V B (A:) ) x 

<5(k-k'), (46) 

(P ± (k)P;(k')>-(27T)- 2 (P £ (fc)-P B (fc))x 

e ±«V>*<5(k-k'). (47) 



(43) 
(44) 
(45) 



Just as in the scalar case, the visibility associated with 
a pointing center b and a baseline u can be expressed in 
terms of the Fourier transform of the antenna pattern: 

V±(u,b) = (2tt) 2 J d 2 fci*(k + 2vru)P ± (k)e 4k - b . (48) 

The correlation between two visibilities is 



V 



±± 

12 



"12 



(2tt) 2 y d 2 fci*(k + 27rui)i(k + 27ru 2 ) 
^■^- b2 \V E {k)+V B {k)), (49) 

{^f J d 2 fci*(k + 27rui)i(k + 27ru 2 ) 
e lk < hl ~ h2 \V E (k) - V B (k))e ±4 ^ k . (50) 



5.2. Spherical harmonics 

Since the quantities P± are quantities of spin weight 
=p2, it is natural to expand them in spin- (±2) spherical 
harmonics: 



Pi(f) = ^ a^2,lro T 2Y em {r). 



(51) 



The decomposition into E and B components is par- 
ticularly simple in terms of the spherical harmonic coef- 
ficients: 

a±2,im = E( m ± iBi m . (52) 

The two-point statistics are completely described by two 
power spectra Cf E , Cf B : 



(Ei m Eg, m ,) —Cf E Su'S. 
(Be m Bg, m ,) — Cf B 5u'8. 

(ElmB^ m ,) =0. 



mm' ; 
ram' j 



(53) 
(54) 
(55) 



As in the case of temperature anisotropy, the spherical 
and flat-sky power spectra are related via Ci ~ P(w) 
with I — 2nu. 

The visibilities can be expressed in terms of the spher- 
ical harmonic coefficients as 



a=F2,£mP=F 2/m(u)i 



where 



F: 



=p2/m 



(u) 



d 2 r A(v) T2 Y im (r)e 2 



(56) 



(57) 



Consider first the covariance between two visibilities 
with identical pointing centers. Combining equations 
(|39|) and (|52|l through H56fl . the visibility covariances can 
be shown to be 



V 



±± 



=E(c? + cf)wf t I 



where 



i 

± = ^P T 2/m(ui)P 



=f2/m 



(ua). 



(58) 
(59) 

(60) 



In the case where the two observations have different 
pointing centers, we once again transform to a coordinate 
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Fig. 6. — Polarization window functions. The left panel shows the window functions for the correlation (VqVq) of the two baselines at 
the upper corners of Fig. [5] The right panel shows the correlations calculated by incorrectly assuming the entire mosaic is flat — that is, 
ignoring the rotation of the basis vectors 0, <j> when moving from one point to another. Note that the vertical axes of the upper plots differ 
by a factor of 200. These window functions were calculated in the cylindrical approximation, but the exact spherical harmonic calculation 
yields negligible differences. 



system with both pointing centers on the equator, sepa- 
rated by an angle {3. Because the spin-weighted spheri- 
cal harmonics have azimuthal dependence e lm ^, the only 
change is an additional factor of e lm ^: 



and 



irn(3 



(61) 



In order to calculate the correlation between a pair of 
observations with arbitrary pointing centers, we simply 
rotate to a new coordinate system that places both cen- 
ters on the equator before applying the above results. 
In performing this rotation, it is important to remem- 
ber to transform P± (and hence V±) by e ±2lS where S 
is the angle through which the polarization basis direc- 
tions are rotated by the transformation. To be specific, 
if the change of coordinates results in a rotation of the 
6, <j> directions at each of the pointing centers by Si, S2, 
then V12 

explicit recipe for finding these angles. 

5.3. Connecting flat-sky to spherical. 

As in the case of temperature anisotropy, we can see 
the connection between the spherical and flat-sky calcu- 
lations of polarization by considering observations that 
lie near the equator of our spherical coordinate system 
and approximating the sphere by a cylinder. By app ly- 
ing the spin-raising operator fe.g.. lLewis et al.ll2002(l to 
the plane- wave approximation to the spherical harmonics 
<|2U|) . one can show that in this limit 



V^e 2 *^* 1 ^). See the Appendix for an 
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By reasoning similar to the previous section we can use 
this to connect the spherical harmonic formalism to the 
flat-sky limit. 

5.4. Example 

Consider a 10 x 10 mosaic of pointings of an interfer- 
ometer, with the same parameters as in the example of 
Section the Gaussian beam width is a = 5°, and the 
pointings are centered on the equator and separated by 
5° in both 9 and </>. We consider only one baseline per 
pointing, with u = 220. Assume that both visibilities 
Vq and Vjj are measured (either directly or by measur- 
ing V±). 

For any pair of pointings we can define E and B win- 
dow functions such that 

( V Qi V Qi) = Y,« C * E + W ? G ? B ) ( 65 ) 

1 

where Vqi is the visibility for Stokes Q corresponding to 
baseline i. In the limit of infinitely sharp ^-space reso- 
lution, we would expect Wf to vanish, since the polar- 
ization would be parallel to the baseline u. Inevitably, 
however, when only part of the sky is covered (leading to 
imperfect Fourier space resolution), there i s some mix- 
ing of E and B modes l|Lewis et al.ll2002l iBunr] 12003 
iBunn et alJEOQ^ . 

For the case of the two visibilities at the upper corners 
of the grid, these window functions are shown in Fig. El 
As in the scalar case, the correlations are strongly af- 
fected by the rotation of the coordinate basis. If we in- 
correctly model the entire mosaic as flat, treating the 
basis vectors 6, <p at each point to be parallel, the cor- 
relation between these two pointings would be dramati- 
cally overestimated. In fact, by treating the sky as flat, 
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Fig. 7. — Window functions for a single pointing (solid), the sum 
of all pointings (dashed), and the sum of all pointings neglecting 
sky curvature (dotted). 



we would be making two separate errors: treating the 
two baseline vectors as parallel (just as in the case of 
temperature anisotropy in Sec. |3J) and failing to apply 
the appropriate transformation to the Stokes parameters 

(Q, u). 

Fig. [7| illustrates the improvement of resolution due 
to mosaicking in this example. Like Fig. [3| this figure 
shows the autocorrelation window function for a single 
pointing as well as that of the sum of all 100 pointings 
in the grid. If we neglect sky curvature, we overestimate 
the correlation between distant baselines and hence also 
overestimate the improvement in £-space resolution. 

6. CONCLUSIONS 

Interferometers have been used to great effect in mea- 
suring CMB temperature and polarization anisotropies. 



The formalism for analyzing interferometer data, how- 
ever, has only been fully developed in the small-field-of- 
view or flat-sky limit. Future experiments which aim for 
exquisite £-space resolution will need to survey large ar- 
eas of sky - outside the realm of validity of the existing 
formalism. 

In this paper we have extended the formalism to the 
situation where we can approximate the sky as flat for 
each individual pointing of the instrument, but we re- 
lax the assumption that the angle between pointings is 
also small. We have connected the full-sky spherical har- 
monic approach to the flat-sky Fourier approach in two 
distinct ways and derived approximations for the visibil- 
ity covariance matrix in each. We find that the cylin- 
drical method of and H5.2I works in all cases better 
than the harmonic method of and provides accurate 
approximations to the full-sky expressions for individual 
pointings smaller than 20° FWHM. Mosaicking together 
many pointings increases the i?-space resolution, but in 
the cases considered here the improvement is less than 
would be predicted from the flat-sky formalism, in large 
part due to the effects of baseline rotation. If we ne- 
glect sky curvature we overestimate the correlation be- 
tween distant baselines and hence also overestimate the 
improvement in f-space resolution. 
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APPENDIX 

In calculating the covariance between visibilities at two different pointing centers ?i,f 2 , we must transform to a 
coordinate system that places both pointing centers on the equator. This affects both the components of the baseline 
vectors Ui,U2 and (in the case of polarization) the Stokes parameters Q,U. We present here an explicit recipe for 
performing this transformation. 

Throughout this appendix, unprimed symbols will refer to the original coordinate system, and primed symbols will 
refer to a coordinate system (x' , y' , z') such that both f i lies on the x' axis and r 2 is in the x'y' plane. First, choose 
the z' axis to be perpendicular to both vectors: 

* = ^ (Al) 
|ri x r 2 | 

Next, choose the y' axis to be perpendicular to both z' and fj.: 

y' = f i x z'. (A2) 
Finally, choose x' = y' x z'. In spherical coordinates (8', <fi') defined with respect to the primed coordinate system, we 
have fx = (n/2,0) and f 2 = (tt/2,/3) with {3 such that (cos (3, sin j3) = (f 2 • x',f 2 • y'). 

Say that the baseline vector (z = 1,2) is expressed in the original (unrotated) spherical coordinate system as 

Ui = u lS 9 + u^cf). (A3) 

We need to know the corresponding components in the rotated coordinate system. The components of the basis vectors 
6, 4> in the rotated coordinate system are 

e = {cos§i)e'- (sin (A4) 

4> = (sin 8i)0'+ (cos 5i)(j>'. (A5) 
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The easiest way to find the rotation angle Si is to compute the components of <j> = (z x ri)/|z x fj|, and take the dot 
product sine),; = cf> ■ 9' = <p • {— z') since is on the equator in the primed coordinate system. 
Once the rotation angles 8\ , 62 are known, the components of the baseline vectors are 



Ui0i \ _ f cos Si sin Si \ / u i0 
Uip ) y-sinJi cos Si J \Ui4, 

In calculating the polarization visibilities we replace (Q,U) with 

Q'\ _ ( cos2(5 4 sm.25i\ ( Q 
U' I { -sin2(5 4 cos2(5 4 ) ' I U 



(A6) 



(A7) 
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